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IN VIVO BRAIN ELASTICITY MEASUREMENT BY MAGNETIC 
RESONANCE ELASTOGRAPHY WITH VIBRATOR COIL 

CROSS-REFERENCES TO RELATED APPLICATIONS 
[0001] NOT APPLICABLE 

STATEMENT AS TO RIGHTS TO INVENTIONS MADE UNDER 
FEDERALLY SPONSORED RESEARCH AND DEVELOPMENT 
[0002] NOT APPLICABLE 

REFERENCE TO A "SEQUENCE LISTING," A TABLE, OR A COMPUTER 
PROGRAM LISTING APPENDIX SUBMITTED ON A COMPACT DISK. 
[0003] NOT APPLICABLE 

[0004] This invention relates to magnetic resonance elastography (MRE). Specifically, a 
15 vibrator coil imparts vibration to the brain during MR. The imparted vibrations allow a non- 
invasive determination of brain tissue elasticity for diagnosis of patient risk for malignant 
brain edema and herniation following acute brain trauma, stroke, intra-cerebral hemorrhage 
and other brain disease. 

BACKGROUND OF THE INVENTION 

20 [0005] The mechanical properties of biological tissue often change during pathological 
processes. This is probably nowhere more evident that in neurological disorders, where as a 
consequence of brain encasement in the rigid skull vault, my brain tissue swelling may have 
little room for compensation. Disturbances of brain tissue elastance or compliance by trauma, 
stroke, infection or neoplasm results in alteration of the intra-cranial pressure and pressure- 

25 volume curve. (Thompson W. Intracranial Hypertension. In: Oh T, ed. Intensive Care 

Manual. London: Butterworths, 1990.) Compensatory treatment may result in a decrease in 
the intra-cranial blood volume, a decrease in CSF volume or osmotic shrinkage of more 
elastic brain tissue using mannitol or hypertonic saline osmotic pressure gradients. 
[0006] Acute stroke affects about 500,000 people a year (Wolf P, D f Agostino, RB. 

30 Epidemiology of Stroke. In: Barnett H, Mohr, JP, Stein, BM, Yatsu, FM., ed. Stroke. New 
York: Churchill Livingstone, 1998:3 - 28.) and is the third most common cause of death in 
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15 



20 



25 



30 



the United States. About 10% of all strokes involve occlusion of the middle cerebral artery 
with hemispheric infarction (Bogousslavsky J, Van Melle, G, Regli, E. The Lussanne Stroke 
Registry: Analysis of 1000 consecutive patients with first stroke. Stroke 1988; 19:1083- 
1085; Sacco R, Toni, D, Mohr, JP. Classification of Ischemic Stroke. In: Bamett H, Mohr, 
JP, Stein, BM, Yatsu, FM., ed. Stroke. New York: Churchill Livingstone, 1998:341- 400.). 
About 10% of these patients develop coma and malignant brain edema (MBE) (Melo T, de 
Mendonca, A, Crespo, M, Carvalho, M , Ferro, JM. An emergency room-based study of 
stroke coma. Cerebrovascular Diseases 1992; 2:93 - 101); with an associated mortality of 
-80% (Bushnell C, Phillip-Bute, BG, Laskowitz, DT, et al. Survival and outcome after 
endotracheal intubation for acute stroke, Neurology 1999; 52:1374 - 1380.). MBE is 
secondary to tissue swelling due to increased cell water content following ischemia and 
cellular metabolic failure. Such swelling may results in herniation of brain tissue by a sub- 
falcine, transtentorial, tonsillar or rostral-caudal mechanism causing compression and 
compound damage to non-ischemic brain tissue. In the case of ischemic damage, maximal 
brain swelling usually occurs within 3-5 days post onset of stroke. Another cause of stroke 
related abnormal brain swelling and edema are intra-cerebral bleeds. 
[0007) Head injury and associated brain trauma are also major public health problems 
being a major cause of mortality and morbidity in the 1-44 year age group (Kraus J, 
McArfhur, DL, Silverman, TA, Jayaraman, M. Epidemiology of Brain Injury. In: Narayan R, 
Wilberger, JE, Povlishock, JT., ed. Neurotrauma. New York: McGraw-Hill, 1996.). After 
head trauma brain swelling may result in unequal brain compartment pressure gradients with 
resultant tissue shifts and herniation. Clinically the patient often demonstrates a deteriorating 
level of consciousness together with more localizing neurological signs. The timing of best 
medical therapy and surgical decompression in brain herniation syndromes is unclear. In 
some categories of stroke such as ischemic cerebellar ischemic stroke . or hematoma surgical 
decompression has clearly been shown to be advantageous (Heros R. Cerebellar hemorrhage 
and infarction,. Stroke 1982; 13:106 - 109; Jauss M, Krieger, D, Horning, C, Schramm, J, 
Busse, O. Surgical and medical management of patients with massive cerebellar infarctions: 
results of the German-Austrian Cerebellar Infarction Study, Journal of Neurology 1999; 
246:257 - 264.) . Other more investigative decompressive techniques such as 
hemicraniectomy and duroplasty have been shown in uncontrolled patient series to reduce 
patient mortality and morbidity following acute hemispheric stroke (Deleshaw J, Broaddus, 
WC, Kassell, NF, Haley, EC, Pendleton, GA, Vollmer, DG, Maggio, WW, Grady, MS. 
Treatment of right hemispheric cerebral infarction by hemicraniectomy. Stroke 1990; 21 :874 
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- 881; Carter B, Oglivy, CS, Candia, GJ, Rosas, HD, Buonanno, F; One-year outcome after 
decompressive surgery for massive nondominant hemispheric infarction. Neurosurgery 
1997;40:1168-1176.). 

[0008] The optimal timing of intervention either medical or surgical is often unclear with 
5 clinical indices often having a low sensitivity and specificity in predicting timing of 

intervention (Schwab S, Steiner, T, Aschofif, A, Schwarz, S. Steiner, HH, Jansen, O, Hacke, 
W. Early hemicraniectomy in patients with complete middle cerebral artery infarction. 
Stroke 1998; 29:1888 - 1893) . 

[0009J MRE is a relatively recently implemented MR technique enabling non-invasive 

10 measurement of tissue elasticity by imaging alteration of the magnetic spin density caused by 
mechanical vibration or displacement wave propagation through deeper tissue with 
amplitudes in the order of a few micrometers. The technique is developing but has previously 
been used to examine breast and muscle tissue. By defining tissue elasticity MRE may 
provide unique imaging information of acute brain herniation syndromes allowing the design 

15 of more applied clinical studies where these questions could be systematically studied 

(Muthupillai R, Lo mas , DJ, Rossman, PJ, Greenleaf, JF, Manduca, A, Ehman, RL. Magnetic 
resonance elastography by direct visualization of propagating acoustic strain waves, Science 
1995; 269:1854 - 1857 ; Muthupillai R, Rossman, PJ, Lomas, DJ, Greenleaf, JF, Riederer, SJ, 
Ehman, RL. Magnetic resonance imaging of transverse acoustic strain waves. Magnetic 

20 Resonance in Medicine 1996; 36:266-274; Van Houten E, Miga, MI, Weaver, JB, Kennedy, 
FE, Paulsen, KD. Three-Dimensional Subzone-Based Reconstruction Algorithm for MR 
Elastography; Magnetic Resonance in Medicine 2001 ; 45:827-837; Weaver J, Van Houten, 
EEW, Miga, MI, Kennedy, FE, Paulsen, KD. Magnetic resonance elastography using 3D 
gradient echo measurements of steady-state motion. Med. Phys. 2001; 28:1620-1628; Van 

25 Houten E, Paulsen, KD, Miga, MI, Kennedy, FE, Weaver, JB. An overlapping subzone - 
technique for MR-based elastic property reconstruction. Magnetic Resonance in Medicine 
1999; 42:779-786). 

BRIEF SUMMARY OF THE INVENTION 
[00101 A vibrator coil is applied to the skull by adaptation of a commercially available 
30 transcranial Doppler monitoring harness during MR applies mechanical waves, in Oie 
auditory acoustic range, through the skull to the brain, typically at the temporal acoustic 
window of the skull. Utilizing magnetic resonance elastography (MRE), non-invasive 
estimation of tissue elastic properties in three dimensions occurs. The propagation of the 
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acoustic waves through brain tissue results in transverse phase alteration of voxel isochromats 
allowing measurement of brain elasticity in the presence of applied magnetic field gradients. 
A protocol of timing for the acoustical excitation of the brain in a range from 125 hertz to 500 
hertz is disclosed which includes synchronizing the acoustical interrogation to the subject's 
5 heart beat with a period of pre-excitation of the brain before the gating the interrogating radio 
frequency to the head of the patient. Image processing of the final data received from the MR 
scan includes image accumulation of phases in opposing directions, subtraction of the 
accumulated images to obtain a phase map, unwrapping of the phase map, to extract absolute 
phase and finally relating the phase to the original displacement to obtain the elastic 

1 0 properties of the brain being examined. 

[001 1] The clinical motivation for such measurements is to determine normal brain tissue 
compliance and the pathological alteration of brain tissue compliance or elasticity. Brain 
tissue compliance or elasticity alterations occur in neurological conditions such as brain 
trauma, acute stroke associated with malignant cytotoxic edema and in brain tumors 

1 5 associated with vasogenic edema. Tissue swelling results in altered mechanical properties 
while continued tissue swelling progressing to brain herniation often results in a reduced 
patient functional outcome. The present procedure leads to the non-invasive measurement of 
both the normal and altered brain compliance enabling the norm to be identified, compared to 
the abnormal, and allow identification and timely intervention for some of the above 

20 neurological conditions. 

BRIEF DESCRIPTION OF THE DRAWINGS 
I0012J Fig. 1 is a perspective view of a patient having the vibration inducing coil mounted 
to his head; 

[0013] Fig. 2 is an enlarged perspective view of the vibration inducing coil; 
25 [0014] Fig ;3A is a plan view of the probe in contact with the acoustic window of a patient's 
head; 

[0015) Fig 3B is a planned view of the spring biased probe from a mounting bracket 
illustrating the vibration inducing coil position for imparting acoustical vibration to the probe; 
[0016] Fig 3C is a side elevation of the probe with the coil not shown illustrating the 
30 preload bar for applying adjustable torque to the probe at the acoustic window of the patients 
head; 

[001 7] Fig 3D is a side elevation of Fig. 3C; 
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[0018] Fig 4 schematically illustrates the synchronous trigger pulses, the motion-sensitizing 
gradient, the interrogating RF pulse, and the three-dimensional acquisition of data here 
schematically illustrating data acquisition of front to back sections taken through the skull; 
[0019]^ Fig 5 is a block diagram of the required processing for signals received from the ■ 
5 MR device for determining the magnetic resonance elastography (MRE) of the brain being 
examined; > 

[0020] Fig. 6 is a block diagram of the sequence of triggering the coil with respect to the 
MRI machine; 

[0021 J Fig. 7 illustrates a patient at the head with a required "birdcage" about the head 
10 illustrating the collar with the vibration inducing coil mounted to his head within the volume 
defined by the bird cage, this illustration illustrating in the background the MR tunnel with 
the vibrator coil of Fig. 3A-3C hidden from view;, and, 

[00221 Fig. 8 A to 8D illustrate a MRE single slice technique on a healthy volunteer 
showing sagittal slice acquisition with 
15 [0023] Fig. 8A illustrating excitation frequency 125 Hz; 

[0024] Fig. 8B illustrating Hilbert Transform; 

[0025] Fig. 8C illustrating phase unwrapping of the Hilbert Transform; and, 

[0026] Fig. 8D illustrating a Shear modulus map. 

DETAILED DESCRIPTION OF THE INVENTION 
20 [0027] Referring to Figs. 1 and 2, a small custom-built MR compatible coil 10 is placed 

about 5 cm from the temporal window so that torque, perpendicular to the static MR Bo static 
field occurs after passing a small alternating current through the coil 1 0. The coil uses the Bo 
static MRI magnetic field and the applied current to cause vibration of a fulcrum 20 against 
the temporal window 31 of skull 30. This results in the generation of displacement waves 
25 within the skull and intra-cranial cavity causing displacement of tissue isochromats. 
Vibration frequencies in the range of 125 - 1000 Hz are used. The coil is comfortably 
applied to the skull 30 by adaptation of a commercially available transcran ial Doppler 
monitoring harness 40. Because of the presence of the ambient magnetic field B 0 , coil 10 
actuates without its own self contained magnet. 
30 [0028] The standing wave field of mechanical stress required for MRI Elastography is 

produced by the vibrations of the small lightweight coil 10 mounted to an available ultrasonic 
transcranial Doppler device made MR compatible into holder 40 that attaches to the head as 
shown in Fig L The centerline of the coil is oriented at right angles to the B 0 static magnetic 
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field, and when an alternating current from a remote, computer controlled power amplifier 
passes through the coil (not shown), it rotates about an axis 1 1 (shown in Figure 2). This axis 
1 1 is perpendicular to both B 0 and the coil centerlinel2 (See Fig. 3A). 
[0029] The coil is rigidly attached to a shaft 12 supporting the coil 10. A rigidly attached 
5 contact rod 14 touches the skull at the acoustic window 31 just in front of the upper ear. The 
coil 10 is designed so that a sinusoidal coil current at, for example, 250 Hz, results in a 250 
Hz vibration applied to the skull. A torsion spring 16 set in a preload mechanism is attached 
to the coil shaft 12 allowing the contact rod 14 to be pushed against the skin with an 
adjustable and operator defined force through preload bar 18 and spring 19 to ensure optimal 
10 skull apposition. An accelerometer (not shown) attached to the contact rod can allow the 
vibration amplitude to be monitored. Several adjustments are available to insure that the 
contact rod is able to touch the skull at the acoustic window for different sized people while 
still maintaining the coil in the proper orientation with respect to B 0 . 

[0030] The coil 10 consists of 80 turns of 1/4 mm diameter copper wire wound in a single 

15 layer onto a 3 mm thick, 5 cm diameter by 2.5 cm long plastic cylinder. An excitation 

current of 2 amperes at a frequency of 500 Hz (most demanding case) is predicted to produce 
an angular oscillatory amplitude of 4.5 milliradians. The corresponding motion at the end of 
the 3 cm long contact rod is approximately 0.1mm when it barely touches the skin, and will 
be somewhat less when a preload is applied. Since this motion is applied at approximately 45 

20 degrees to the surface, both compressive and shear stresses are applied to the skin. The 
accelerometer output divided by the excitation frequency provides a signal proportional to 
velocity and can be used to determine how to vary the current at different frequencies. 
[0031J In practice, the mounting bracket 42 will be adjusted in the head frame 41 so that 
the contact rod 14 lines up with the acoustic window 31 in the anterior posterior direction and 

25 tightened. Nextthe head frame 40 will be attached and tightened; With the coil aligned - 
perpendicular to B 0i contact rod 14 is rotated about the coil shaft 12 until it touches the skin at 
the acoustic window 31 and it is tightened on the shaft. (This may require adjustment of the 
coil support together with further tightening,) Finally the preload will be adjusted to the 
desired value through preload bar 18 and spring 19. 

30 [0032] Previous work (by Muthupillai R, Lomas, DJ, Rossman, PJ, Greenleaf, JF, 
Manduca, A, Ehman, RL. Magnetic resonance elastography by direct visualization of 
propagating acoustic strain waves. Science 1995; 269:1854 - 1857 and Muthupillai R, 
Rossman, PJ, Lomas, DJ, Greenleaf, JF, Riederer, S J, Ehman, RL. Magnetic resonance 
imaging of transverse acoustic strain waves. Magnetic Resonance in Medicine 1996; 36:266- 
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274) has demonstrated the principle of MRE in agarose phantoms and in human muscle. We 
adopt these techniques allowing direct visualization of acoustic strain waves by synchronized 
MR imaging with motion sensitizing gradients. Development of the NMR signal acquisition 
code takes place on a 1.5T Signa GE machine. A standard gradient echo sequence will be 
modified so that synchronized vibration-coil and gradient echo images are obtained (See Fig. 
3 of the referenced article). 

[0033J By using transmissible acoustic waves through the cranial cavity MR phase changes 
are obtained. Since they are sound waves in the range of 125 - 1000 Hz, they have no 
adverse biological consequences. A point from a sinusoidal signal generator will be used to 
trigger coil vibrations to produce a mechanical steady state before application of the initial 
MR radio frequency pulse and motion-sensitizing gradients. The initial phase of the 
symmetric motion-sensitizing gradients will vary by ±x The following image acquisition 
parameters are typical: pulse repetition time (TR) cardia gated -1000 ms for a subject with a 
heart rate of 60 beats/minute ms, echo delay time (TE) 24 ms, slice thickness 5.0 mm, 128 - 
256 phase encoding views with an image acquisition time of about ~ 90 s, NEX -2-4, 
gradient field of- 900 - 3500 mGauss. Duty cycle of the coil 10 will be - 5% of the TR time. 
This is calculated by multiplying the motion-sensitizing period T by the number of cycle and 
dividing by TR. 

[0034] Fig. 7 illustrates patient observation during acoustical excitation of skull 30. This 
figure illustrates patient with the apparatus of Fig. 1 attached placed within standard MR 
"birdcage" 50 and passed into MR apparatus 60. Fig. 8 is representative scans brain scans of 
volunteer individuals taken utilizing the process of this invention. 

[0035] Referring to Fig. 5, software written in National Instruments Lab VIEW is used to 
output the excitation waveform that drives the excitation coil 70. The software creates the 
waveform and controls a digital-to-analog converter 72 (National Instruments PCI-6070E) 
that outputs the waveform. A current stabilized amplifier 74 is used to generate the currents . 
necessary to drive the excitation coil. A highpass filter 73 is used to prevent DC currents 
from driving the coil (see Fig. 5). 

[0036] The software allows the user to excite the coil before the actual MRI acquisition 
begins. The waveform was designed to minimize currents during the period of the MRI 
sequence where the RF pulse is output and during data acquisition. It was found that any 
currents present during these parts of the sequence produced artifacts in the MRI images. 
[0037J Data Analysis 
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[0038] Refeiring to Fig. 5, a block diagram of the data analysis here present is illustrated. 
MR P files 90 are processed by fast Fourier transform (inversion) software. Phase 
unwrapping then occurs to eliminate ambiguities of phase signal redundancy. As will be set 
forth below, three discrete types of processing are possible utilizing Hilbert transforms 95, 
local wave length 96 and shear modulus 97. What foUows is a theoretical explanation of 
these computer confined techniques. 

10039] A magnetic-field gradient results in a phase shift $ of the NMR signal and is given 
by 

l 004 °] #-rJo r *X/)* 

o 

[0041] where y is the gyromagnetic ratio for a proton, Gr is the magnetic gradient field, t is 
the time duration of the gradients while r(t) describes the position of the nuclear spins as a 
function of time. The following analysis assumes that the tissue properties are locally 
isotropic with the stress-strain relationship defined by a Hookean relationship. A stationary 
vibration field will be allowed to evolve prior to sampling the magnetization field by the 
receiver coil. Considering the case where the acoustic source is periodic and applied for a 
sufficient period to damp transients then the deviation from equilibrium of an isochromat is 
determined by local wave properties. These are modified from the source by attenuation, 
reflection and scattering resulting in a mixture of transverse and longitudinal waves 
determining the local dyadic strain tensor values 

[0042] The local stress-strain relation is given by the following tensor equation 30 
[0043] F=A,(y £ )l + 2fiE 

since the time varying body force Fl is 

[0044] - ^.V-F_=/> 0 ~. . 

[0045] which becomes 

25 [0046] (^)VV* +/ ,V.V*= A |£ 

[0047] by considering the irrotational field the above equation becomes 

[0048] cf a* 

[0049] and by consideration of the solenoidal field the wave equation becomes 
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[0050] 
[0051] 

[0052] 
[0053] 

[0054] 
[0055] 

[0056] 
[0057] 
ratio. 
[0058] 



V.Ve= 



1 &s 



where 



, L l -2a- J 



and 



and 



2(1+<t) Po 



c2 ^ Y(l-a) 
' A (l+o-)(l-2cr)/7 0 

hence it possible to isolate for each voxel the Young's modulus and the Poisson 



The spin density p(r) is a vector field and related to the NMR signal in the receiver 
coil. The variation in the spin density is proportional to the strain tensor in a manner 
dependent on the local stress tensor generated by the acoustic field where 

s„ 0 0 
0 



0 



p(r)cc 

[0059] 

following the assumption of local, isotropy. 

[0060] By subtracting the spin density p with the applied motion sensitizing gradient 
acoustic strain field from the spin density with the motion sensitizing gradient opposite in 
phase, the phase representation of the strain field can be derived. The above analysis makes 
the assumption of tissue isotropy. this can be further extended by taking the trace of the 
strain tensor resulting in a mean local strain sav. 



Tr 



20 [0061] 



0 



yy 
0 



0 
0 



[0062] This can be further simplified to give the average strain for the whole brain. 
Discrimination of white and gray matter (WM, GM) might be possible using a histogram 
fitting technique so that average WM and GM properties could be defined. The magnitude of 
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MR gradient is set to be synchronized with the acoustic strain field induced by the HeadWave 
coil and may described by 

[00631 = ^ G ' ;/ e {"T&n+yr/zoriin+iyrnjin+iyr} 

[0064J where n=0, 1,2. . .N-l, and T=2tt/g>. The phase shift in the received signal from the 
stationary displacement field induced in the tissue by the coil vibration is given by 

[00651 * 1 ^ 

[00661 where k is the wavenumber.r is the distance from the source, a is the phase lag 
between the MR gradient and the vibration coil and £ is the local tissue displacement. A 
displacement field can be obtained after inverse Fourier transformation of the k-space 
complex image and subtraction of the two out of phase motion sensitized gradients images. 
The reconstructed phase subtraction image allows derivation of the wavelength X and 
wavenumber (k) (A.=2n/|kD by taking the Hilbert transformation (See Fig. 5 at 95) and 
calculating the instantaneous phase followed calculation of the shear modulus G. The 
Young's modulus may be calculated from the shear modulus and by assumption of tissue 
isotropy (See Fig. 5 at 97). By taking images in the two other orthogonal directions, the tissue 
properties can be isotropically determined. Further images taken in the cross directions 
allows anisotropic determination of tissue properties. 

[00671 The equation describing the displacement in an isotropic medium is the Navier 
equation as follows 

[0068] k mj u p ^ 

[00691 where X and u are the Lame constants. This can be solved directly as a forward 
problem. The inversion of the 3D isotropicjwaveequation solving for X and u, can be.solyed. 
using a least squares minimization method 15, 17. 

[00701 A tensor, Eeff of Young's moduli (E) may be derived with MR by altering the 
gradient combination (X, Y, Z, X+Y, X+Z, Y+Z). Eeff is a symmetric second order tensor 
with the dominant directions and magnitudes determinable from the tensor eigenvalues and 
eigenvectors on a pixel-by-pixel basis. The solution of the full 3D wave equation for the 
velocity vector field followed by calculation of the full Young's tensor Eeff adds a further 
order of complexity not only in the data analysis but also in the data acquisition where the 
complete tensor strain field must be obtained. 
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[0071] Outcome Measures 

[0072] Primary Outcome Measure 

10073] The phase difference images (<f>) represent the priniary outcome measure. Increasing 
levels of complexity will be developed from single slice single plane data to 3D volume 
acquisition in three orthogonal directions. The most complex level of data acquisition 
possible consists of tensor acquisition in six independent directions in each plane. 
100741 Secondary Outcome Measures 

[0075] The following data will be derived from the phase difference images 
[0076] [1] Wavelength (X) of the acoustic wave in different areas of the tissue together with 
derivation of the shear modulus (n) where n=(v^) 2 p (v is the excitation frequency, p is the 
medium density) 

[0077] [2] Displacement field (§) 

[0078] [3] Derived material properties in terms of the Lame coefficients for an isotropic 
medium 

[0079] Statistical Analysis 

[0080] Determination of sample size (power analysis) 

[0081] We are determining some indication of the random error involved in MRE 
following simulation experiments using a 2% agarose phantom skull 30 when multiple 
measurements are made. This random error most likely will be ah underestimate of the 
random error obtained across a patient population since the variance from measurement in 
patients will be added to the variance inherent in MRE measurement. Further the biological 
variance is likely to be the dominant component. 

[0082] We are undertaking a healthy volunteer, investigation to acquire a normal population 
data. At the present time, the total number of subjects required is unclear until data analysis 
gives some indication of the data variance invol ved in MRE. By the central limit theorem, 
analysis of population sizes of -30 or greater should approach a Gaussian distribution 
allowing the application of normal statistics to MRE data analysis. This will be of importance 
in allowing standard brain imaging analysis techniques such as SPM (statistical parametric 
mapping). SPM fully takes into account multiple comparisons procedures by the application 
of Gaussian fields. The normal volunteer data acquired will allow appropriate power 
calculations to be generated in the evaluation of subsequent clinical hypotheses. - 
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[0083] Methods used to analyze outcomes 

[0084] The following statistical considerations will be applied to analysis of the primary 
and secondary outcome variables 

[0085] [1] Each subject will be treated as independent so that single global subject 
5 dependent variables will be compared by normal univariate statistics. 

[0086] [2] All multiple variables from single subjects will have appropriate multiple 
comparison techniques applied. For imaging data, this will entail that each of the independent 
scalar component of either, the derived vector or tensor fields from the primary or secondary 
outcome variables will be analyzed using Gaussian field theory and statistical parametric 
10 methods. Statistical parametric methods take into account necessary multiple comparison 

corrections and the data will be analyzed with SPM, a standard statistical parametric package. 
[0087] Validation of MRE in Healthy Volunteers 

[0088] Following the above data analysis the values derived from MRE can be compared 
with the experimentally determined Lame coefficients and Young's moduli for central 

15 nervous tissue available in the literature (See Hagermann A, Rohr, K, Stiehl, HS, Spetzger, 
U, Gilsbach, JM. Biomechanical modeling of the human head for physically-based, non rigid 
image registration. IEEE TRANS MED IMAGING 1999; 1 8:875 - 884.). 
[00891 Fig, 8 A through 8D are graphic depictions of a MRE single slice technique on 

a healthy volunteer showing ipsilateral sagittal acquisition. The MR imaging was performed 

20 using a cardiac gated, phase contrast, gradient echo sequencel .5T, TE 26 ms, FOV 18 cm, 
256 x 128, slice thickness 5.0 mm with motion-encoding gradients (7 to 2 cycles, 3.5 G/cm) 
applied during the TE period. All phase images are windowed with land marking slice 
acquisition through the point of actuator apposition in the axial plane. In sagittal slice 
acquisition the motion encoding gradients are in the frequency direction. 

25 [0090] Fig. 8 A illustrates an excitation frequency 125 Hz, Gradient 3 .5G/cm, windowed to 
±1 .4 radians and phase unwrapped showing areas of high and low phase accumulation 
secondary to the transmitted transverse acoustic wave, pre-excitation 117 msec. This is a 
schematic of the phase measurement in a sagittal brain slice after phase unwrapping and 
windowing to +/- 1.4 radians. 

30 [0091] Fig. 8B illustrates a Hilbert Transform of Fig. 8A showing the 90° phase shift. This 
allows the mathematical generation of the 90 degree quadrature image. 
[0092] Fig. 8C illustrates phase unwrapping of the Hilbert Transform of Fig.8B. After 
further phase unwrapping the Hilbert transform and combining Fib. 8A and 8B, the 
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instantaneous phase can be derived on a pixel-by-pixel basis. By differentiating the local 
phase we get the local frequency or wavenumber from which the shear modulus can be 
calculated as in Fig. 8D. 

[0093] Fig. 8D illustrates the shear modulus map derived from the local spatial frequency 
or wavenumber map calculated by differentiating the instantaneous phase and applying 
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WHAT TS CLAIMED IS: 

1 . A method for magnetic resonance elastography of at least a section of 
the brain comprising the steps of: 

examining the head of a patient in vivo in a magnetic resonance device; 
vibrating the head of the patient during the examination at a selected 
frequency between 125 hertz and 500 hertz; 

observing and plotting phase alteration of voxel isochromats at the selected 
frequency to obtain phase patterns; and, 

measuring the phase patterns across at least the section of the brain. 

2. The method for magnetic resonance elastography of at least a section 
of the brain according to claim 1 wherein the measuring by observing phase patterns 
includes: 

repeating the examining, vibrating, observing and plotting, and measuring 
steps for a group of individuals; and, 

comparing the measuring of the phase patterns from one individual to other 

individuals. 

3. The method for magnetic resonance elastography of at least a section 
of the brain according to claim 1 wherein the measuring by observing phase patterns 
includes: 

analyzing the phase patterns utilizing Hilbert transforms. 

4. The method for magnetic resonance elastography of at least a section 
.of the brain according to claim 1 wherein the measuring by observing phase patterns 
includes: 

analyzing me phase patterns by utilizing the shear modulus. 

5. The method for magnetic resonance elastography of at least a section 
.of the brain according to claim 1 wherein the measuring by observing phase patterns 
includes: 

analyzing the phase patterns by utilizing the local wavelength. 
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6. The method for magnetic resonance elastography of at least a section 
of the brain according to claim 1 wherein: 

the observing and plotting phase alteration of voxel isochromats occurs after 
vibrating the head of the patient for about a time period of 5 - 200 msec. 

7. A method for magnetic resonance elastography of at least a section of 
the brain comprising the steps of: 

affixing a coil to the head of the patient in a magnetic resonance device having 
a magnetic field; 

passing alternating current through the coil to cause vibrational energy to pass 
from the coil to the head of the patient at a selected frequency between 125 hertz and 500 
hertz; 

after the passing step, examining the head of a patient in the magnetic 
resonance device; 

observing and plotting phase alteration of voxel isochromats at the selected 
frequency to obtain phase patterns; and, 

measuring the elasticity of the brain by observing the phase patterns across at 
least the section of the brain. 

8. The method for magnetic resonance elastography according to claim 7 

and wherein: 

observing and plotting phase alteration of voxel isochromats at the selected 
frequency to obtain phase patterns immediately after passing of the alternating cuirent 
through the coil has ceased but before vibrational energy within the head of the patient 
dissipates. . " ~ " 

9. The method for magnetic resonance elastography according to claim 7 
and wherein the affixing of a coil to the head of the patient includes: 

placing a shaft through the coil to receive vibrations from the coil; 

placing a probe in contact with a shaft at one portion and biasing the probe 
into contact with a human skull at another portion; and, 

vibrating the coil to impart vibrations through the shaft to the probe to vibrate 
in vivo a human brain within the skull. 
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10. The method for magnetic resonance elastography according to claim 9 
and wherein the biasing of the probe into contact with the human skull includes: 

biasing the probe into contact with the oral cavity of the human skull. 

11. A method for magnetic resonance elastography of at least a section of 
the brain comprising the steps of: 

examining the head of a patient in vivo in a magnetic resonance device; 

observing the periodicity of the patient's heartbeat for determining a sampling 
interval with respect to the patient's heartbeat; 

vibrating the head of the patient immediately before a sampling interval at a 
selected frequency between 125 hertz and 500 hertz; 

observing and plotting phase alternation of voxel isochromats at the selected 
frequency to obtain phase patterns; and, 

measuring by observing the phase alternation across at least the section of the 

brain. 

12. The method for magnetic resonance elastography of at least a section 
of the brain according to claim 1 1 comprising the further steps of: 

ceasing the vibrating immediately before the observing and plotting step. 

13. An apparatus for improved magnetic resonance analysis of the brain 
during magnetic resonance examination comprising: 

a mounting for biasing a probe on to the cranium of the patient in a magnetic 
resonance device; 

a coil affixed to the probe for passing vibrations from the coil to the probe; 

and, 

means for passing an alternating current through the coil in the range of 125 
hertz to 500 hertz to cause the coil to vibrate within the magnetic field of the magnetic 
resonance device and pass the vibrations of the coil to the probe. 
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14. The apparatus for improved magnetic resonance analysis of the brain 
during magnetic resonance examination according to claim 13 and wherein means for passing 
alternating current through the coil includes: 

a high pass filter and a current stabilized amplifier. 
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Patent 

Attorney Docket No.: 015280-484000US 



IN VIVO BRAIN ELASTICITY MEASUREMENT BY MAGNETIC 
RESONANCE ELASTOGRAPHY WITH VIBRATOR COIL 

ABSTRACT OF THE DISCLOSURE 
[0094] A vibrator coil is applied to the skull by adaptation of a commercially available 
transcranial Doppler monitoring harness during MR applies mechanical waves in the acoustic 
waves through the skull to the brain. Utilizing magnetic resonance elastography (MRE), non- 
invasive estimation of tissue elastic properties in three dimensions occurs. The propagation 
of the acoustic waves through brain tissue, coupled to phase alteration of voxel isochromats 
in the presence of applied motion encoding magnetic field gradients allows measurements of 
brain elasticity. 

SF 1428849 vl 
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